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We consider tlie asymptotically flat double-Kerr solution for two equal mass black holes with 
either the same or opposite angular momentum and with a massless strut between them. For fixed 
angular momentum and mass, the angular velocity of two co-rotating Kerr black holes decreases 
as they approach one another, from the Kerr value at infinite separation to the value of a single 
Kerr black hole with twice the mass and the angular momentum at the horizons merging limit. We 
show that the ratio J /KP for extremal co-rotating Kerr black holes varies from unity at infinite 
separation to two at the merging limit. These results are interpreted in terms of rotational dragging 
and compared with the case of counter-rotating Kerr black holes. We then analyse the merging of 
ergo-regions. In the co-rotating case the merger point occurs at an angle of 7r/2, in agreement with 
recent general arguments. In the counter-rotating case the ergo-regions never merge. We study the 
horizon geometry for both cases as a function of the distance and provide embedding diagrams. 
Finally, we study the thermodynamical evolution of the co-rotating double Kerr system, showing 
that, in the canonical ensemble, it is thermodynamically stable for fast spinning black holes. As 
for single Kerr black holes the stable and unstable phases are separated by a second order phase 
transition. We show that for large fixed angular momentum two Kerr black holes reach a minimum 
distance, before horizon merging has occurred, where the thermodynamical approximation breaks 
down. We also consider the micro-canonical ensemble to study the maximal energy that can be 
extracted from the double Kerr system as a function of the separation between the black holes. 



I. INTRODUCTION 



Understanding the interaction between two Kerr black 
holes is of crucial importance in connecting the General 
Theory of Relativity with observations of astrophysical 
black hole binaries . One expects that the full dynamical 
evolution of a binary system of Kerr black holes, in par- 
ticular the merger stage, can only be tackled in numerical 
relativity. Despite the remarkable recent progress, which 
has allowed one to infer the properties of the final state in 
binary black hole coalescence (see [l| for a review), there 
are still issues which the current numerical codes have 
not yet addressed. For instance, the codes do not keep 
track of the horizon angular velocities of individual black 
holes, but only of their angular momentum 2]. The an- 
gular velocities are of physical interest because they are 
associated to the existence of an ergo-region which, in 
turn, is associated to potentially observable phenomena 
such as superradiance. 

As a first step for understanding this and other phe- 
nomena in a realistic binary system, one can consider a 
quasi-static approximation, where relative velocity effects 
are neglected and the black holes are kept in equilibrium 
by an external force. Studies of such systems include 
the double-Schwarzschild solution 0, 0, [a] and the dou- 
ble Kerr solution 6]. The latter is an exact solution of 



•Electronic address: 'miguelc@fc.up.pt 
^Electronic address: crherdei@fc.up.pt 
■f Electronic address: imrebelo@fc.up.pt 



the four dimensional vacuum Einstein equations which, 
for some region in parameter space, describes two Kerr 
black holes rotating in the same plane, in an asymptot- 
ically flat spacetime and with a massless strut between 
them, which provides the equilibrium. For instance, we 
have recently shown (3] that for a special case of the dou- 
ble Kerr system (for equal Komar masses Mi = M = M2 
and opposite Komar angular momenta Ji = J = — J2) 
one can provide an exact formula in terms of the physical 
parameters for the angular velocity of the black holes: 
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where ( is the coordinate distance between the black 
holes in Weyl canonical coordinates. This angular ve- 
locity varies from the Kerr value as ^ -f-oo, to zero 
as C ^ 2M, corresponding to the touching of the two 
horizons. Such slow down of the angular velocity is very 
intuitively interpreted as a consequence of the rotational 
dragging that the black holes exert on each other. 

As a consequence of ([T]), one may suspect that, in this 
counter-rotating system, the ergo-regions never merge 
because, as one decreases the distance between the black 
holes, they will become closer and closer to the horizons, 
collapsing onto the horizons when the horizons them- 
selves touch. This shall be explicitly shown in this pa- 
per. One should note that in this system the Kerr upper 
bound for the angular momentum depends on C 
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varying from the Kerr value as ^ +00, to infinity 
as C ^ 2M. Thus, the ergo-regions never touch even 
tough the individual black holes can have an arbitrar- 
ily large angular momentum, for fixed mass, if they are 
sufficiently close. The "trick" is that, despite the arbi- 
trarily large angular momentum, the angular velocity is 
always sufficiently small as to allow a black hole hori- 
zon. This example, therefore, makes crystal clear that 
it is the angular velocity, and not the angular momen- 
tum, that determines the existence of a horizon, for an 
uncharged, rotating black hole, an intuitive observation 
if one invokes the membrane paradigm 

If one considers the co-rotating system (equal Ko- 
mar masses Mi = M ~ AI2 and angular momenta 
Ji — J — J2), one might naively expect an opposite 
behaviour. That is, as they approach one another, the 
mutual rotational dragging should speed up the horizon 
angular velocity of either black hole. Reasoning along the 
same lines as in the last paragraph, one could then expect 
that the extremality limit would decrease as compared to 
the usual Kerr value. These expectations, however, im- 
mediately lead us into trouble. A first argument is the fol- 
lowing. If this were true, two under-extreme Kerr black 
holes (no strut between them), dropped from a certain 
distance with zero velocity could become extreme and 
over-extreme as they would fall into one another. Indeed, 
the initial distance and angular momenta could be tuned 
such that velocity effects are as small as one wishes, so 
that one could trust the results of the double-Kerr solu- 
tion as a quasi-static approximation. This would clearly 
violate cosmic censorship and, potentially, the second law 
of thermodynamics. A second argument is that, in this 
co-rotating case, one expects a single Kerr black hole with 
twice the mass and twice the angular momentum to arise 
when the black holes merge. This will indeed be true, the 
reason being clearly understood from the rod structure 
of the co-rotating double Kerr solution, to be described 
below. For a single Kerr black hole, the angular velocity 
as a function of mass and angular momentum obeys 
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Thus, if the double Kerr co-rotating solution interpolates 
between two Kerr black holes of mass M and spin J and 
one Kerr black hole of mass 2M and spin 2 J, the angu- 
lar velocity of the horizon must slow down as we decrease 
the distance between them. Consequently, the extremal- 
ity limit should increase as compared to the usual Kerr 
value. It is simple to see that a co-rotating double Kerr 
system that interpolates (fixing M, J) between a final ex- 
tremal black hole and two initial Kerr objects, requires 
the latter to obey J = 2M^. Thus, we expect the ex- 
tremality limit to vary between the Kerr value, asymp- 
totically, and J = 2Af^ just before merging. 

The analysis of the co-rotating double-Kerr system is 
more involved and no simple expressions analogue to ([T]) 
and (21) could be obtained. Nevertheless, we will exhibit 
all physical quantities in a way analogous to the ones 





FIG. 1: Horizon embedding for the counter-rotating (left) 
and co-rotating (right) double Kerr system, for two different 
values of tlie physical distance d and for M = 1 = J. 



of the counter-rotating case. We shall show numerically 
that, for fixed mass and angular momentum, the angular 
velocity of the two black holes indeed slows down as they 
approach; and that the ratio J/M^ for extremal black 
holes grows from unity, asymptotically, to two, in the 
merging limit. It follows that, in this co-rotating case, 
we expect the ergo-regions of the two black holes to in- 
deed merge, and eventually become a single one, as they 
approach. 

The behaviour obtained in the co-rotating system can 
still be interpreted in terms of rotational dragging. The 
point is that each black hole is simultaneously dragging 
and being dragged. To drag some material object effec- 
tively increases the moment of inertia of the black hole 
and hence decreases its angular velocity. In the counter- 
rotating Kerr system, the dragging exerted by and upon 
each black hole work with the same purpose: slowing 
down each one of them. But in the co-rotating sys- 
tem they work with opposite goals: on the one hand, 
dragging the other black hole slows down the black hole 
which is dragging; on the other hand, being dragged by 
the other black hole speeds up the black hole which is 
being dragged. The observed behaviour shows that it is 
the first effect that is the dominant one, and this was 
exactly the effect that was being neglect in the naive ex- 
pectations mentioned above. This interpretation is qual- 
itatively confirmed by the fact that, for the same M, J 
and physical distance d, the angular velocity is smaller in 
the counter-rotating, where both effects slow down each 
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FIG. 2: Horizon angular velocity in terms of the physical dis- 
tance for M — 1 and J = 1.25, 1.5, 5/3, for both the counter- 
rotating (grey) and co-rotating (black) cases. We have chosen 
values of J/M^ greater than unity to enhance the difference 
between the two cases. For these choices the solutions de- 
scribe black holes only up to some maximal physical distance. 



horizon, than in the co-rotating case - Fig. [51 

Since in the co-rotating case we expect the two ergo- 
spheres to merge we may ask what is the merger angle 
at which this happens. Recently, a universal behaviour 
for mergers of ergoregions has been observed 0. For 
D + 1 dimensional vacuum solutions of the Weyl class 
with rotation in a single plane, the merger angle is always 
given by 



9m = 2arccot-\/(5 — 1 , 

where S — D — p and the merger is extended along p 
spatial dimensions. In our case D = 3. To determine 
the value of p we recall Hajicek's theorem [l^l for ergo- 
surfaces in four dimensional General Relativity: at the 
fixed points of the C/(l) rotational isometry, ergo-surfaces 
must either touch the horizon or hit the black hole singu- 
larity. Thus, for the double Kerr solution, we expect that 
at the symmetry axis {p — in Weyl canonical coordi- 
nates) the ergo-spheres will only touch when the horizons 
touch. Hence the merger should be at p > and there- 
fore p — 1 due to the U{1) symmetry. Thus, the expected 
angle at the merger point is 6'm — 7''/2 for the double Kerr 
system, and we shall indeed show this to be the case. 

In the co-rotating Kerr solution one can use standard 
instanton methods to study the thermodynamics of the 
system since, unlike the counter-rotating case, both the 
temperatures and the chemical potentials (angular veloc- 
ities) of the two black holes match. For the single Kerr 
solution it is well known [11] that there is a second order 
phase transition at Jc/M'^ — 0.68 or, in terms of the tem- 
perature and angular velocity, at ilc/T = 5.85. At this 
phase transition, the specific heat at constant J has an 
infinite discontinuity, being negative for slowly rotating 
black holes (f2 < fie) and positive for fast spinning ones 
(fi > ric)- We shall show that in the double Kerr system 
there is a similar behaviour. In particular, when con- 
sidering the canonical ensemble, fixing the total angular 



momentum, the double Kerr system is thermodynami- 
cally stable for fast spinning black holes. The system is 
then described by its temperature, total angular momen- 
tum and distance between the black holes. We shall see 
that, provided the angular momentum is large enough, 
the black holes will evolve to approximate one another, 
until a minimum distance is reached. At this critical 
point, which occurs before the horizons merge, the spe- 
cific heat diverges and the thermodynamical approxima- 
tion breaks down. One then expects large entropy, and 
therefore horizon area, fluctuations to occur. In this situ- 
ation it is a very interesting question to understand what 
happens to the system, even before the horizons merge. 

Another question that can be analysed in the frame- 
work of gravitational thermodynamics is that of the max- 
imal energy that can be extracted from the collision of 
two Kerr black holes, which occurs for adiabatic pro- 
cesses. In the case of co-rotating black holes, we shall 
confirm the bound of 29% derived by Hawing [Tj] , study- 
ing the dependence of this energy on the angular momen- 
tum and distance between the black holes. 

This paper is organised as follows. In section |TT] we 
shall describe the rod structure of the Kerr solution and 
its ergo-sphere in Weyl canonical coordinates, since these 
will be used for the double-Kerr solution. In particu- 
lar, we shall develop an intuition about the shape of the 
ergo-sphere in these coordinates. We shall moreover de- 
scribe how one can construct embedding diagrams for 
the horizon. In section Hill we briefly describe the double 
Kerr solution and its rod structure. In section IIII Al we 
analyse the counter-rotating case. We shall review the 
main results of 0] and introduce a suggestive notation. 
Moreover, we show that the ergo-spheres of the two in- 
dividual counter-rotating Kerr black holes never merge 
and we construct the embeddings of the horizon in 3- 
dimensional Euclidean space as shown in the diagrams 
of Fig. [1] These embedding diagrams give an intuition 
about the horizon deformations due to the strut and the 
interaction between the two black holes. In section HlI Bl 
we analyse the co-rotating Kerr black holes case in a 
form which is similar to that of the counter-rotating case. 
Moreover, we exhibit the angle dm = tt/2 at the ergore- 
gions merger point. We plot numerically the behaviour 
of the angular velocity of the horizon of both black holes 
as a function ( for fixed M, J, for various values of J, 
which exhibits the "slow down" effect we have described 
above. A numerical plot of the extremality condition 
Jext/M^, as function of the coordinate distance C is also 
given, again exhibiting the behaviour discussed: for fixed 
M, Jext/M"^ increases from unity, as we decrease C,, to 
two. In section IIVI we analyse the thermodynamics of 
the co-rotating double Kerr system. We start by com- 
puting the action of the double Kerr instanton, therefore 
computing the free energy in the grand-canonical ensem- 
ble. We then consider the canonical ensemble and anal- 
yse thermodynamical stability. Finally, we move to the 
micro-canonical ensemble to study energy extraction in 
the process of merging two Kerr black holes. 
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FIG. 3: Rod structure for the single Kerr solution. 

II. KERR SOLUTION 

In this paper we shall analyse four dimensional, asymp- 
totically flat spacetimes in Weyl canonical coordinates 
(i, (f), p, z). The line element is of the form 



gttdt^ + 2gt4,dtd(i) + g^^dcj)"^ + e^^^P^'-Xdp"^ + dz^) 



where gab — gab{p,z). The (f) coordinate has periodicity 
2tt and is chosen such that gt^, — > at infinity. That 
is, at infinity we have a non-rotating frame. Thus, ergo- 
surfaces, if they exist, will always be defined by gu = 0. 

Weyl solutions are conveniently described by their rod 
structure, including the direction of the rods The 
rod structure of the Kerr solution is given in Fig. [3l The 
solution may be obtained by a 2-soliton transformation, 
at z = ai (anti-soliton) and z — a2 (soliton) with BZ 
vectors 

(l,2ai6) and (l,2a2c) , 

respectively, from the Schwarzschild solution. The an- 
tisymmetric choice b = — c immediately guarantees that 
the NUT charge is zero and therefore that the solution is 
asymptotically flat. The result is a two parameter family 
of vacuum solutions, with parameters given by the length 
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It turns out that in terms of the physical parameters 
of the solution 
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and M and J = Ma denote respectively the ADM mass 
and angular momentum of the Kerr black hole. 

For future reference it is also useful to recall that the 
area and temperature of the Kerr black hole are com- 
pactly written as: 
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Note that AT = 021. This is a generic result that also 
holds for individual black holes in the double Kerr solu- 
tion. A similar statement is true for Smarr's formula 
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FIG. 4: Kerr black hole ergo-surface in Weyl canonical coor- 
dinates. We fixed M = 1 and varied the angular momentum 
J = aM. 



It follows from Smarr's formula that the angular velocity 
of the outer Kerr event horizon is 



2Mr+ 

In Weyl canonical coordinates the Kerr ergo-sphere is 
defined by the equation 
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where 
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In Fig. |3]we show the location of the Kerr ergo-sphere in 
Weyl canonical coordinates, for various values of a. The 
location of the (spatial sections) of the outer horizon is, 
on the other hand, given by p = and ai < z < a2- 
Introducing a polar coordinate 9 by 

2z = oi + 02 + 021 cos 9 , 

the induced metric on the horizon becomes 
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The conformal factor is constant, A ~ 2mr+, and 
/ ^ n 2n 1 + ^^ 

1 + O^M^ 

The above 2-geometry can be embedded in flat Euclidean 
3-dimensional space using the embedding functions 
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FIG. 5: Profile of the Kerr black hole horizon embedding, 
for fixed mass M — 1 and different values of the angular 
momentum J — aM. 



as long as g'{u)'^ < 4 (prime denotes derivative). This 
always holds in the equatorial plane, for any Kerr black 
hole; but does not hold for sufficiently fast spinning Kerr 
black holes and for sufficiently high latitudes In Fig. 
Owe exhibit the profile of the embedding, for fixed mass 
and different values of the angular momentum. 



III. DOUBLE KERR SOLUTION 

The double Kerr solution Q is a seven parameter 
family of solutions of the vacuum Einstein equations. 
Physically the seven parameters may be written as 
Ml, M2, Ji , J2, Ci Maxis and bNUT, which are respectively 
the two black hole masses and angular momenta, the co- 
ordinate distance between the two black holes, the mass 
of the axis between them and the NUT charge (see [7J 
for a more detailed description). Here we shall focus on 
asymptotically flat solutions {hMUT — 0), obeying the 
axis condition {Maxis = 0) and with equal mass for the 
black holes M = Mi = M2 ■ Moreover we shall either con- 
sider J = Ji — J2, the co-rotating case or J = Ji = — J2, 
the counter-rotating case. This leaves us with two 3- 
parameter families of solutions (M, J, () . For both these 
families, the rod structure is displayed in Fig. [6l The 
coordinate distance is defined as 



C = ^32 



043 + 021 
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We shall refer to the black holes with horizon at ai < 
z < a2 and 03 < z < 04 as the first and second black 
holes, respectively. 
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FIG. 6: Rod structure for the asymptotically flat double Kerr 
solution, obeying the axis condition, for equal mass black 
holes with equal (fii = ^2) or opposite (f2i = — ^2) spin. 



A. Counter-rotating Double Kerr 

This case was first considered in [l^ and studied in de- 
tail in 0] . The solution is obtained by a 4-soliton trans- 
formation, at 2: = fli (anti-soliton), z = 02 (soliton), 
z = 03 (anti-soliton) and z = 04 (soliton), respectively 
with BZ vectors 

(l,2ai6) , (l,2a2c) , (l,2a3c) and (1,2046) , (8) 

from the double-Schwarzschild solution. Moreover, we 
choose from the beginning 

021 — 043 • 

This symmetric choice of BZ vectors and rod configura- 
tion immediately guarantees that the axis condition is 
obeyed. Additionally, asymptotic flatness is guaranteed 
if one takes 



(c + 6)(l + 5c) . 
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The result is a three parameter family of vacuum solu- 
tions, with parameters given by the BZ parameters h and 
c and the coordinate distance Q (or equivalently 032 , due 
to (0)). 

For a physical analysis it is desirable to take as inde- 
pendent parameters the Komar mass M of either black 
hole, the Komar angular momentum J (in modulus) of 
either black hole and the physical distance d between the 
black holes. Using the physical distance complicates mat- 
ters way too much to sustain any hope of finding simple 
analytical expressions. So, we shall contempt ourselves 
with taking as independent "physical" variables M ,J and 
the coordinate distance Q. This is justifiable since C is a 
good measure of the physical distance 



d(C,M, J) 



in the sense that, for fixed M and J, d grows monoton- 
ically with C, as shown in Fig. [71 Therefore we use the 
latter, since it simplifies the analysis. 

In analogy to the single Kerr solution, we introduce 
the constants 



M± v/M2 - a2/(C, M) 



where M denotes the Komar mass of either black hole, 
J = Ma is (minus) the Komar angular momentum of the 
first (second) Kerr black hole and 



/(CM) 



C- 2M 
C + 2M 



6 



d 

50 r 



{M,J)=(2,2) (M,J)={1,1) 

/ / / / 

{M,J)={2,4)/ // /{M J)={l]0.5) 



10 - 




FIG. 7: Relation between coordinate distance C, and physical 
distance d in the counter-rotating double Kerr solution, for 
fixed values of M, J . The (dashed) straight line d = is given 
for comparison. For J jM^ < 1, the slope of d = d(^) always 
approaches unity, asymptotically. 



It turns out that in terms of these quantities, the BZ 
parameters can be simply re-expressed as 



a/ (CM) 



the constraint ^ becomes 

and ([7]) takes the form 

a32 = C - 2^AP-a^f{C,M) . 
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It follows from Smarr's formula and the above results 
that the angular velocity of the outer Kerr event horizon 
is given by 



2Mr- 



-/(C,M) = -fi2 , 



(12) 



which is relation ([T]) , given in the introduction. Since the 
extremal limit is obtained as 021 — 043 — > 0, it follows 
from pn)) that it is given by 



ext 



which is relation ([2]), also given in the introduction. 

Finally, the area and temperature of the two black 
holes can be compactly expressed as 
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FIG. 8: Force between the two black holes for two physical 
distances d = 3 and d = 5 as a function of J for Af = 1, for 
the counter-rotating (grey lower curve for each distance) and 
co-rotating cases (black upper curve for each distance). 
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Notice that the introduction of the constants r± displays 
a remarkable similarity between the physical quantities in 
this counter-rotating double Kerr solution and the ones of 
the single Kerr solution exhibited in the previous section. 

It follows from pTjl that ( 2M is the horizons touch- 
ing limit, since 032 0. Naively it seems that in this 
limit the rod structure in Fig. [6] reduces to that in Fig. [3] 
with fl+ = 0, by virtue of (fT2|l . In other words, that the 
double Kerr solution with fixed M and J interpolates be- 
tween two counter-rotating Kerr black holes (at C 00) 
and a single Schwarzschild black hole (at C = 2AI). In- 
deed, in the latter limit, the temperature (fT4|) yields that 
of a single Schwarzschild black hole with mass 2M, and 
the total area of the two black holes Ai + A2 (fT3)) yields 
that of a single black hole with mass 2M. However, in 
general, this is not the case, because the BZ vectors of 
the solitons at z = 02 = 03 do not become trivial, cf. 
([8]). Thus, only if one takes J — > 0, the touching limit 
becomes a merging limit, and one recovers the (single) 
Schwarzschild solution; otherwise the resulting geometry 
is singular. 

The force between the two black holes, computed from 
the massless strut between them, is exactly given by 
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Naively this is independent of J. However, if expressed in 
terms of the physical distance d, rather than the coordi- 
nate distance the force will depend on J. One expects 
that the spin-spin force should be attractive if the black 
holes are counter-rotating (see Q for a discussion of this 
point). Thus, one expects the modulus of the force to 
increase as a function of J, for fixed physical distance. 
Although this is a naive expectation (cf. discussion in 
the next section) it is indeed verified in Fig. [51 
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FIG. 9: Double Kerr ergo-surface in Weyl canonical coordi- 
nates. We fixed M = 1, J = 0.8 and varied the coordinate 
distance as indicated. 

In Weyl canonical coordinates, the location of the ergo- 
surface is determined by (taking into account ([9])) 

{{r+~r^fx+Hf + l)]'rj}'-[a{,f + l){r+-r^)e]' = 0, 

(16) 

where x, rj and e are the following functions of p and z: 

X = Ml3Ai24P?2P34 1 

V = Ml2Ai34P?3P24 , (17) 

e = hm^sP [m24Pi3 - Aiis^y > 
and we have used the shorthand notation 

/Xy = - , pfj = + PtPj . (18) 

In Figs. [51 and [TUl we display the location of the counter- 
rotating double Kerr ergo-surface. From the plots it is 
clear that the ergo-spheres of the two individual black 
holes never touch. 

As for the single Kerr case one can show that the in- 
duced metric on the horizon of each black hole takes the 
form ([5]). One can then use functions (O to construct 
an embedding for the horizon. In Fig. [11] we display the 
profile of such embeddings for the first black hole and for 
various physical distances (fixing M and J). One notices 
that, unlike the single Kerr case, the embedding is not 
symmetric, as expected from the existence of a strut in 
one of the sides (upper), as well as the interaction with 
the other black hole. Since we have chosen, M = 1 = J, 
these are extremal black holes at infinity. Thus, the em- 
bedding covers only a small part of the horizon. However, 
as they approach, they slow down considerably and the 
size of the embedding surface increases. In particular, in 
the "south pole" , where there is no strut, the horizon be- 
comes completely covered by the embedding. By making 
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FIG. 10; Double Kerr ergo-surface in Weyl canonical coordi- 
nates. We fixed M — 1 and varied both the angular momen- 
tum and the coordinate distance (as indicated) such that as 
to obey the extremality condition 



a revolution of these profiles we obtain the 3-dimensional 
embeddings displayed in Fig. [TJ 




FIG. 11: Profile of the horizon embedding for the first 
("lower") black hole in the counter- rotating (grey lines) and 
co-rotating (black lines) double Kerr systems. We take 
M — J = 1 and various values for the physical distance 
d = 75,25, 10,5,3, 1. 
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FIG. 12: / as a function of C for M = 1 and J = 0.8, 1.2, 1.5. 
The detail shows a zoom near the merging points ^o- The 
higher the angular momentum, the lower the maximum al- 
lowed value for 

B. Co-rotating Double Kerr 
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FIG. 13: Angular velocity, for AI — 1 and fixed J in terms of 
the coordinate distance in the co-rotating double Kerr solu- 
tion. 



Two equal Kerr-Newman black holes were discussed 
in [T^ , where the solution was generated via transforma- 
tions of the Ernst potential and the issue of force balance 
was addressed. Herein we generate it using the inverse 
scattering construction and focus on different physical 
issues. 

The solution is obtained by a 4-soliton transformation, 
at z = ai (anti-soliton) , z — a2 (soliton), z = (anti- 
soliton) and z = 04 (soliton), respectively with BZ vec- 
tors 

(l,2ai6) . (l,-2a2c) , (l,2a3c) and (1,-2046) , 

(19) 

from the double-Schwarzschild solution. Moreover, we 
choose from the beginning 

021 = 043 ■ 

This antisymmetric choice of BZ vectors on this symmet- 
ric rod structure immediately guarantees that the NUT 
charge is zero and therefore that the solution is asymp- 
totically flat. Additionally, the axis condition is obeyed 
provided 

(l + 2c2 + 6V)^ + 2cif^^ = i-l + bcf . (20) 
C b~c C 

The result is, again, a three parameter family of vacuum 
solutions, with parameters given by the BZ parameters b 
and c and the coordinate distance C (or equivalently 032 , 
due to ([7])). Note that the relation between coordinate 
and physical distances is qualitatively very similar to the 
one for the counter-rotating case displayed in Fig. [71 

In this case we were not able to express b and c in 
terms of the "physical" parameters (M, J, Q . The rea- 
son can be traced back to the fact that (PO)) gives an ex- 
pression for 021 = 021(6, c, C) less tractable than the one 



obtained from ^ . Replacing back in the expressions for 
the Komar quantities one obtains M — M{b,c,(^) and 
J = J(6, c, ^), which could be inverted in the counter but 
not in the co-rotating case. 

Nevertheless, a similar structure to that of the counter- 
rotating case can be displayed, in terms of an analogue 
function /(^, A/, J), which we determine numerically. 
That is, defining 

where M and J = Ma denote the Komar mass and Ko- 
mar angular momentum of either black hole, the phys- 
ical quantities are expressed by identical expressions to 
the ones in the last section. Thus, the constraint pp)) 
becomes 

021 = a43 = f+ - = 2^M^-a^f{C,M, J) , (21) 
and ([7]) takes the form 

032 = C - 2^Jm^ -a^fiCM, J) . (22) 

The merging limit corresponds to 

C ^ Co = V(2M)2-o2 . 

Thus, the behaviour of /(C, M, J) in this limit and in the 
asymptotic limit is 

hm /(C, M, J) = 1 , Hm /((, Af, J) ^ ] ■ 

The interpolating behaviour is displayed in Fig. 1121 
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FIG. 14: Jext/M'^ in terms of the coordinate distance, for 
the co-rotating case and the counter-rotating case. We fixed 
M = 1 which determines the merging limit in the counter- 
rotating case to be at C = 2. 



It follows from Smarr's formula and the above results 
that the angular velocity of the outer Kerr event horizon 
is given by 



2Mr 



^/(C,M, J) =02 



(23) 



The angular velocity is shown graphically in Fig. 1131 
One notes that, for fixed M, the angular velocity O de- 
creases with decreasing distance, as anticipated in the 
introduction. Since the extremal limit is obtained as 
0-21 = 043 ^ 0, it follows from (|21[) that it is given by 



The comparison between the extremal limit in the 
counter and co-rotating cases is given in Fig. [TH One 
observes that for extremal black holes the ratio Jext/M"^ 
increases with decreasing C for both cases, as a conse- 
quence of the rotational dragging slow down. However, 
whereas it diverges in the counter-rotating case, it con- 
verges to two in the co-rotating case, corresponding to 
the Jext /Af ^ value of the Kerr black hole with mass 2M 
and angular momentum 2 J, that is formed after merging. 

The area and temperature of the two black holes can 
be compactly expressed as 



Ai = A2 = 87rA/f+/i(C, M, J) 



SvrAff 



(24) 



(25) 



+ 



The function h is exhibited in Fig. [Tsl Note that it 
has the same asymptotic behaviour, at infinity and in 
the merging limit, as the corresponding function in the 
counter-rotating case, cf. (|13p and 



FIG. 15: ft as a function of C for M = 1 and J = 0.8, 1.2, 1.5. 
The detail shows a zoom near C, = C,q. The higher the angular 
momentum, the lower the maximum allowed value for C,. 



The force between the two black holes, for fixed mass 
and physical distance, is displayed in Fig. [H In this case 
the naive expectation is that the force should decrease in 
modulus, since by general arguments the co-rotating spin 
spin force should be repulsive. Fig. [5] shows, however, 
that this is not the case. The force increases, in modulus, 
as we increase J, for fixed mass and physical distance. 

The breakdown of our naive expectation is, most likely, 
associated with finite size effects. As we vary J, for fixed 
d and M , we are varying the horizon shape and therefore 
higher gravitational multiple moments should contribute 
to the force at the same order as the spin-spin interaction. 
The plot in Fig. [5] is informing us that they mask the 
spin-spin force completely. 

It is, however, reassuring that the force is greater, in 
modulus, in the counter-rotating case than in the co- 
rotating one. Indeed, it has been argued by Hawking 
that more energy can be extracted in a head on col- 
lision of Kerr black holes, if the spins are anti-aligned, 
rather than aligned, with the direction of separation. 
Taking one of the black holes to be much smaller than 
the other. Hawking derived a spin-spin interaction en- 
ergy which was later shown by Wald using a test 
particle (with spin) approximation, to be the work of the 
spin-spin force between the black hole and the spinning 
test particle. Such force is attractive (repulsive) for anti- 
aligned (aligned) spins. 

The fact that our analysis unveils a force which is more 
attractive in the counter-rotating than the co-rotating 
case is in qualitative agreement with the Hawking- Wald 
spin-spin interaction energy. But note that, in any case, 
there is no overlap between their analysis and ours. We 
consider two equal mass and equal spin (in modulus) 
black holes whereas their calculation is limited to a test 
object in the background of a Kerr black hole. This 
means that the finite size effects we are including are 
completely neglected in their approximation. 

In Weyl canonical coordinates, the ergo-surfaces for the 
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FIG. 16: Evolution of the co-rotating double Kerr ergo- 
surface in Weyl canonical coordinates. We fixed M = 1 = J, 
and considered various distances. The ergosurface merging 
happens for C, = C,o = 2.243464, while the horizon merging for 
C, — \/3. The detail in the second figure shows a zoom of the 
merger, in a box of equal length in p and z so that the 7r/2 
angle becomes clear. 

co-rotating double Kerr are defined, in terms of 5, c and 
C, by the following equation (taking into account 

-[bil - c^)y - c{l - b^)x]^ ^ , 

where we are using the definitions p^ . the shorthand 
notation p8)) and 

a = Ail3M24P?3P24 [M23Ail4P^ + P23P14] > 
X = PM23P?4 [M14P22P33 -X~V] , 

y = pijupIs [M23P11P44 -X-v] ■ 

In Fig. [16] we display the merging of the ergo-surfaces 
for (asymptotically) extremal black holes and exhibit the 
merger angle of tt/2. 

As for the counter rotating double Kerr one can show 
that the induced metric on the horizon of each black hole 
takes the form One can then use functions jS]) to 
construct an embedding for the horizon. In Fig. [11] we 
display the profile of such embeddings for the first black 
hole. They are quite similar to the counter-rotating case, 
but the embedding covers a slightly larger portion of the 
horizon, for small distances. The 3-dimensional embed- 
dings are given in Fig. [T] 



IV. THERMODYNAMICS FOR THE 
CO-ROTATING DOUBLE KERR 

For the counter rotating double Kerr the two black 
holes have the same temperature but not the same chem- 
ical potential (angular velocity). Thus, equilibrium ther- 
modynamics can only be meaningfully considered in the 
co-rotating case. Let us therefore study the co-rotating 
double Kerr black hole as a thermodynamical system. 



A. Grand canonical ensemble 

First we consider the grand canonical ensemble with 
temperature T and chemical potential given by the angu- 
lar velocity fl. The free energy can be computed from the 
standard Euclidean quantum gravity approach. One con- 
siders geometries satisfying the periodic boundary condi- 
tion {T,z,p,(j)) = {t + z, p,(j) + if^^), where r is the 
Euclidean time and /3 is the inverse of the temperature, 
and then computes the action of the double Kerr instan- 
ton. The free energy W = W{T, fl) is then related to the 
action by / = PW. 

To find the double Kerr instanton, first one imposes 
the existence of two horizons separated by a coordinate 
length L along the rotational symmetry axis of the Weyl 
canonical coordinates, for instance by choosing the rod 
structure introduced in section IIIII Then, after solving 
Einstein's equations, one finds that the geometry has a 
strut singularity between both black holes. Thus, the 
thermodynamical analysis is performed at fixed coordi- 
nate length L of the strut between both Kerr black holes. 
This variable is related to C and the parameters in the 
rod structure of section IIIII by L = 032 = C ~ 043 ■ It 
should be remarked that L is a mechanical parameter, 
not a thermodynamical one. Indeed, a simple computa- 
tion shows that the entropy density associated to a single 
strut (or to a cosmic string) is zero. The role of this strut 
is simply to impose the mechanical equilibrium between 
both holes (one should not think of thermal fluctuations 
of L, as it is the case with T and fi). With this word of 
caution we shall write bellow W — W{T,n, L), keeping 
in mind that T and Q are thermodynamical variables, 
while L is a mechanical parameter. 

The contribution to the action from both black holes 
arises, as usual, from the boundary term of the Euclidean 
action at infinity and is given by 



There is, however, another contribution to the action 
coming from the conical singularity between the black 
holes, which were kept at fixed coordinate distance L, 
given by 

Ion J on 
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FIG. 17: Free energy W as a function of L for fixed particular 
values of T and Q,. Asymptotically the free energy matches 
the Newtonian gravitational potential energy, while for short 
distances it is finite. 



where 5 = 27r(l — e~'^) is the deficit angle and A = f] e'L 
is the area of the surface spanned by this singularity, 
which is located at p = 0, \z\ < L/2 and < t < /3. The 
function 7 appears in the definition of the metric given 
in and is computed at the strut [25] . Notice that this 
contribution to the action is positive, since the deficit 
angle is negative, as the strut exerts an "outwards" force, 
therefore preventing the black holes from falling into each 
other. Putting everything together we have 

W{T,il,L) = M + ^{1 - e^)L . 

In this expression one should regard both the mass M 
and the function 7 as functions of T, Q and L. For in- 
stance, in the simpler case of the double Schwarzschild 
black hole the dependence of the mass on the tempera- 
ture and coordinate distance is given by the relation Q 



M 



1 L + 2M 



8ttT L + AM 



In a reversible thcrmodynamical process the variation of 
the free energy W is given by 



dW 



where S is the total entropy and is the total angular 
momentum of the system (not the angular momentum J 
of each black hole, as defined in the previous sections). 
As we shall see bellow both quantities are not simply the 
sum of the entropy and angular momentum of each black 
hole, since the interaction between both holes must be 
taken into account. 

Let us now consider a mechanical process where the 
black holes start from an infinite distance and then slowly 
reduce their distance, while their temperature T and 
chemical potential f2 are kept fixed. The mass of each 



black hole Moo, when they are infinitely far apart, is 
given by the single Kerr expression 



1 



The asymptotic value of the free energy Woo — -^00 is 
just the free energy of two isolated Kerr black holes. The 
angular momentum of each black hole at infinite separa- 
tion Joo is also given by the single Kerr value. As usual, 
the force acting on the black holes can be computed from 
the variation of the free energy W — W{T, il, L) with re- 
spect to L 



T =- 



dW 



As the black holes slowly move in, the variation of the 
free energy is given by the work done on the system 



W-Wr. 



TdL' . 



In Fig. [T7]the generic form of the free energy as a function 
of L is shown. The attractive force between the black 
holes is doing work during the process, so that the system 
reduces its free energy W . This free energy is nothing but 
the static interaction potential between the black holes, 
including the spin-spin interaction and the effect of the 
finite size of the black hole horizons (i.e the multipolar 
force terms). Notice, however, that during this process 
the system is exchanging heat with the thermal bath (in 
fact losing heat since the entropy decreases), as well as 
angular momentum, so that the mechanical force here 
computed is not the same mechanical force considered in 
sections IIII Al and IIII Bl For large distances 



W^Mo 



Ml 



correctly reproducing the Newtonian interaction between 
two point-like masses. As expected, the merging of the 
black holes is thermodynamically favoured. The final 
state is then that of a single black hole of mass Moo and 
angular momentum Joo, and the free energy is reduced 
ioW = Moo/2. 

It is interesting to consider the particular case of zero 
temperature, which corresponds to the merging of two 
extremal black holes with fixed angular velocity f2. In 
this case we start with two black holes, each with mass 
and angular momentum given by = Joo = l/(4il^). 
The final state is that of a single black hole with the same 
mass and angular momentum. Thus, just before merging, 
each black hole will have mass Moo/ 2 and angular mo- 
mentum Joo/2. These black holes are extremal and sat- 
isfy the L ~> limiting relation (Joo/2)/(Moo/2)2 = 2. 
Thus, from the thermodynamical evolution at T = and 
fixed f2 we recover the result discussed in the introduction 
and shown in section UlI Bl in fact this is a special case 
of the result therein. Indeed, as we approach two black 
holes varying their Komar mass and angular momentum 
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FIG. 19: Phase diagram (left) and isotherm J — J(f2) (right) 
for the single Kerr black hole. The critical angular velocity 
flc/T — 5.85 corresponds to a critical angular momentum 
obeying Jc/M = 0.68. At this critical line the specific heat 
Cj has an infinite discontinuity, while er changes sign. 



FIG. 18: The entropy as a function of L for fixed particular 
values of T and il. Shown is also the difference between the 
entropy and the contribution from each black hole given by 
the area law. 



as to keep them extremal at all distances, we may, al- 
though we do not have to, keep their angular velocity 
constant. 

The total entropy and angular momentum of the sys- 
tem can also be computed from the thermodynamical 
potential at fixed L, 



S 



dW 
'df 



J 



dW 

'on 



In particular, notice that the entropy computed from this 
expression is not given by the sum of the entropy of each 
black hole obtained from the area law. It contains ad- 
ditional terms, which are subleading at large distances, 
that arise from the interaction between the black holes, 
since the interaction term in W also depends on the tem- 
perature. When the black holes merge we recover again 
the area law. In Fig. [15] we show the entropy of the 
system as a function of the distance L, as well as the dif- 
ference between the entropy and the contribution from 
each black hole obtained from the area law. A similar 
behaviour is observed for the total angular momentum 
JT" of the system, so that only when the black holes are 
infinitely far apart, or when they merge, is the angular 
momentum computed from the Komar integral. 



B. Review of single Kerr thermodynamical 
stability 

To analyse the thermodynamical stability of the dou- 
ble Kerr system, let us first recall some results regarding 
the single Kerr black hole [ll[. In the grand canonical 
ensemble, the stability condition is the positivity of the 



Weinhold metric llSl: its inverse is 



ab 

9w 



dx°-dx^ 



where = (T, VL) (see jl9| for a clear and analogous dis- 
cussion of the thermodynamical stability of the Reissncr- 
Nordstrom black hole). The eigenvalues of this matrix 
are given by /3Cj = /3Co — rf /^t and tr, where 



r -T 



dJ_ 



are, respectively, the specific heat at constant J (which, 
for a single black hole, equals J^) and the isothermal 'mo- 
ment of inertia'. Fig. [T9l illustrates the behavior of both 
Cj and er for the single Kerr black hole: there is a crit- 
ical value of the angular velocity 



T 



5.85. 



below which Cj < and above which < 0, showing 
that in the grand canonical ensemble the Kerr black hole 
is always unstable. 

Working instead in the canonical ensemble, where the 
angular momentum J is fixed, the Kerr black hole has 
a stable phase. In this case one only needs to consider 
the sign of the specific heat Cj, concluding that slowly 
rotating Kerr black holes with < il^ are thermody- 
namically unstable, while fast spinning Kerr black holes 
with > are stable. Both phases, schematically rep- 
resented in Fig. I19[ are divided by a second order phase 
transition, since at the critical value of the horizon an- 
gular velocity the free energy has a discontinuous second 
derivative. This critical point is characterised by an in- 
finite discontinuity of the specific heat Cj. The distinc- 
tion between the two phases can be clearly seen from the 
isotherm J = J(l^), also shown in Fig. [121 Indeed, the 
divergence in the specific heat occurs when the derivative 
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FIG. 20: Isotherms J = J{ft,L) for fixed temperature T 
and three different values of L. The detail on the right shows 
the line of maxima of the isotherms with respect to variations 
of n. This line defines the critical angular velocity Qc — 
n,(T,L). 

of the isotherm J ~ J(f2) vanishes, which in turn imphes 
that the difference between the specific heats at fixed il 
and J diverges. Since Cq is finite at the transition, Cj 
will diverge. 

C. Canonical ensemble 

Clearly, just as the single Kerr black hole, the co- 
rotating double Kerr system is thermodynamically un- 
stable in the grand canonical ensemble. Let us, therefore, 
consider the thermodynamical stability in the canonical 
ensemble, where the angular momentum is held fixed. 
A good starting point is to analyse the behaviour of the 
isotherms = J(yi,L), shown in Fig. [201 We saw in 
section lTV Al that when two black holes start infinitely far 
apart, each with angular momentum J^o = J{^,oo)/2, 
and then move in at fixed T and f2, the final state 
is that of a single black hole with angular momentum 
Joo = Joo(^^,0). It follows that J{n,oo) = 2J^(O,0), as 
can be verified in Fig. [201 Like the single Kerr black hole, 
the double Kerr system also has two distinct phases, one 
of which is stable. To see this, notice that for any fixed 
distance L the isotherm = L) has a maximum 

J^c at some VLc- Hence, fixing the temperature T, as L 
varies from to oo, there is a critical line given by 

with the asymptotic behaviour uj(T, oo) = ^^{T, 0) = 5.85 
given by the single Kerr result. As a function of L we 
observe an oscillatory behaviour in f2c, as shown in the 
detail of Fig. [20l Changing T does not change this qual- 
itative behaviour. We conclude that the specific heat Cj 
has an infinite discontinuity at the critical line, which 
separates the stable, Cj > 0, and unstable, Cj < 0, 



FIG. 21: Behaviour of the specific heat Cj, for a generic 
value of L, as a function of SI. In this plot we set T = 0.1 and 
L = 1, so that the discontinuity corresponds to the maxima 
of the L = 1 curve of Fig. [20l 

phases. As for the single Kerr, both phases are divided 
by a second order phase transition. In Fig. [51] we show 
the behaviour of the specific heat Cj computed for the 
double Kerr system for a finite distance L, which is sim- 
ilar to that of the single Kerr. 

For a given value of the angular momentum J < Jc 
the single Kerr black hole system has two possible states, 
one of which is stable (see Fig. [19]). For the double 
Kerr system this structure is richer since there is also 
the distance parameter L. From Fig. [201 it is clear that 
there are two qualitatively different cases determined by 
the value of the angular momentum J with respect to its 
critical value Jc — Jc{L,T) at L = 0. When the angular 
momentum of the system J satisfies J < Jc{0,T), the 
coordinate distance L may take values from oo to 0. On 
the other hand, when J > J7c(0, T), it is clear that the 
black holes may only reach a minimum distance Lmin ■ In 
both cases there are stable and unstable phases. These 
facts can be seen from Fig. [22l where the slices J = 
const, of the isotherms of Fig. [201 are represented in the 

— L plane. In particular, given J such J > J7c(0, T), 
the black holes will move to the critical line at a finite 
distance Lmm defined by 

To better understand this statement, let us consider the 
thermodynamical potential associated to the mechanical 
evolution of the system at fixed T and J, which is given 

by 

w' = w + nj. 

The black holes will evolve mechanically to minimize the 
potential W = W'{T,J,L) with respect to L. This 
evolution is represented by the arrows in Fig. 1221 In 
Fig. [2Slwe show the behaviour of the potential W' as a 
function of the distance L for the trajectories (a) and (6) 
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FIG. 22: The fixed T and J trajectories in the fl ~ L plane. 
The arrows indicate the mechanical evolution of the system 
given by the thermodynamical potential W' . The red trajec- 
tories (a) describe the merging of fast rotating (top) and slow 
rotating (bottom) Kerr black holes. For the grey trajectories 
(b) both fast and slow rotating Kerr black holes evolve to a 
minimum distance Lmin- Recall that the fast rotating black 
holes are the stable ones. The dashed trajectories correspond 
to the limiting case J = Jc{0,T). 

of Fig. m Clearly, the force 



dL J, J 

is always attractive. 

We have seen that for fixed temperature and fixed and 
high enough angular momentum Kerr black holes have a 
minimum allowed distance\ In the particular case of the 
stable fast rotating black holes, one reasonable question 
to ask is: what happens to the system as the minimal 
distance Lmin is reached? At this point the specific heat 
Cj ^ +00, so that the response of the system to ther- 
mal fluctuations is very large and one can not trust the 
thermodynamical approximation. Fluctuations of the en- 
tropy, and therefore of the horizons area, will be very 
large. For this reason it is tempting to conjecture that 
at this finite distance a phase transition will occur under 
which a single horizon is formed. 

D. Micro-canonical ensemble and energy extraction 

Finally one may also wish to consider the evolution 
of the double Kerr system at fixed entropy S, angular 
momentum and distance L. The corresponding ther- 
modynamical potential is the total energy of the system 
E given by 

E = W' + 

with differential 

dE = TdS + fid J . 
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FIG. 23: Thermodynamical potential W along the trajecto- 
ries (a) and (b) of Fig. [H 
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FIG. 24: Energy as a function of the physical distance d for 
fixed entropy and angular momentum. For both curves the 
asymptotic value of the black holes masses is the same. 



This evolution is relevant to understand the maximal en- 
ergy that can be extracted from the merging of both black 
holes for fixed J', since in that case the variation of the 
energy is maximised by adiabatic processes. In Fig. [24] 
we show the energy of the double Kerr system for fixed 
entropy S and angular momentum J', as a function of 
the physical distance d. Each curve represents a dou- 
ble Kerr system with the same total energy (asymptotic 
masses), but different values of the entropy and angular 
momentum. In both cases we obtain the upper bound 
of 29% for energy extraction derived by Hawking p^ . 
which is independent of the angular momentum. Fig. [24] 
also shows that, for fixed initial energy and until a finite 
distance is reached, one may extract more energy from 
extremal black holes. This fact agrees with the mechan- 
ical analysis of the force between co-rotating Kerr black 
holes given in section ilTTl and illustrated by Fig. [5] since 
the higher J the larger is the force in modulus for fixed 
mass black holes. 
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V. FINAL REMARKS 

In this paper we have studied physical aspects of two 
special cases of the double Kerr solution, which we have 
dubbed the counter-rotating and the co-rotating cases. 
These cases are asymptotically flat, obey the axis condi- 
tion and have Mi = M2 and |Ji| = IJ2I. They simplify 
dramatically the general double Kerr solution, albeit re- 
taining still a considerable degree of complexity. This 
is why we have chosen not to present their explicit con- 
struction herein. The interested reader is referred to 
for details. 

An application of this system is to study rotational 
dragging effects at the level of an exact solution of four 
dimensional General Relativity, as described at length 
in the introduction. A striking physical consequence is 
the possibility that the dimensionless extremality ratio 
Jext/M^ may exceed unity, already discussed in for 
the counter-rotating case. A somewhat analogous phe- 
nomenon takes place for a double Reissner-Nordstrom 
solution of two black holes with opposite charges [1^ [2l| , 
for which the extremality ratio \Qext\/M may exceed 
unity [l^]. A tempting interpretation is that (as for the 
double Kerr solution) it is the electric potential on the 
horizon (angular velocity) rather than the total charge 
(angular momentum) that determines the existence of a 
horizon. It would be interesting to pursue the compar- 
ison between the charged and rotating di-holes in more 



detail. A clear difference is that charged diholes with the 
same charge (see e.g. [2^) yield, in the extremal limit, a 
Majumdar Papapetrou solution, for which the ratio be- 
tween charge (in modulus) and mass is one. 

Another possible application of these double Kerr sys- 
tems is to study geodesic motion in the presence of two 
black holes. This can certainly be done on the equa- 
torial plane of the co-rotating case which (unlike the 
counter-rotating case) provides a totally geodesic sub- 
manifold. A complete study of orbits, however, would 
require the geodesic problem to be Liouville integrable. 
This might sound unlikely. But it was recently pointed 
out by Will [2^ that the Newtonian version of the dou- 
ble Schwarzschild system does posses a hidden symmetry. 
It would therefore be of interest to see if an irreducible 
Killing tensor exists in both the double Schwarzschild 
solution and the particular cases of the double Kerr we 
have considered. 
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